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We examine the extraction of moments of parton distributions from lattice data, focusing 
in particular on the chiral extrapolation as a function of the quark mass. Inclusion of 
the correct chiral behavior of the spin-averaged isovector distribution resolves a long- 
standing discrepancy between the lattice moments and experiment. We extract the in- 
dependence of the valence u — d distribution from the lowest few lattice moments, and 
discuss the implications for the quark mass dependence of meson masses lying on the 
p Regge trajectory. The role of chiral symmetry in spin-dependent distributions, and in 
particular the lattice axial vector charge, g^i, is also highlighted. 



1. Introduction 

Parton distribution functions (PDFs) parameterize fundamental information on the 
nonperturbative structure of the nucleon, and the workings of QCD at low energy. 
Over the past two decades considerable experience has been gained with studies 
of PDFs within low energy models of the nucleon. Ultimately, however, one would 
like a more rigorous connection of PDFs with QCD, and presently this can only be 
provided through the lattice formulation of QCD. 

Because PDFs are defined as light cone correlation functions, it is not possible 
to calculate them directly on the lattice in Euclidean space. Instead, one calcu- 
lates matrix elements of local twist-two operators, which are related through the 
operator product expansion to moments of the PDFs. For the spin-averaged quark 
distributions, q{x) = q'^ (x) + q^{x), the moments are defined as: 

{x"), = f dx a:" {q[x) + {-ir+'q{x)) , (1) 

while moments of the helicity distributions, Aq{x) = q^ {x) — q'^{x), are given by: 

(Ax"), = r dx x" {Aq{x) + (-l)"Aq(x)) . (2) 
Jo 

A number of lattice calculations of PDF moments have been performed over the last 
decade, initially in the quenched approximation, and more recently with dynamical 
quarks. The results indicate that at the relatively large quark masses at which the 
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calculations were made (between w 30 and 190 MeV), the unquenched results 
are indistinguishable from the quenched within the current errors. 

Despite the impressive progress of the lattice calculations, for many years the 
moments have yielded results which were typically ~ 50% larger than the 

experimental values, when linearly extrapolated to the physical quark masses. This 
discrepancy was recently resolved with the observationEI that a linear extrapolation 
in quark mass omits crucial physics associated with the nucleon's pion cloud, and 
that inclusion of the nonanalytic dependence on the quark mass is essential if one 
is to reconcile the lattice data with experiment. 

In this paper we report recent progress made in connecting moments of parton 
distributions calculated on the lattice with experiment, focusing on the extraction 
of the u — d jnoments (which are insensitive to the poorly known disconnected 
contributions^) . We discuss the extraction of the x dependence from moments, and 
the implications for masses of mesons lying on the p Regge trajectory, and outline 
the role of chiral symmetry in helicity distributions and the axial vector charge, gA- 

2. Chiral Extrapolation of Lattice Moments 

The spontaneous breaking of the chiral SU(2)l x SU(2)i^ symmetry of QCD gen- 
erates the nearly massless Goldstone bosons (pions), whose importance in hadron 
structure is well documented. At small pion masses, PDF moments can be sys- 
tematically expanded in a series in mjr, with the expansion coefficients generally 
free parameters. One of the unique consequences of pion loops, however, is the ap- 
pearance of terms nonanalytic in the quark mass, cx m^, which arise from the 
infrared behavior of the chiral loops, and are therefore model independent. 

The leading order (in ttItt) nonanalytic term in the expansion of the moments 
of PDFs was shown by Thomas et alE to have the generic behavior logTOTr- On 
the other hand, in the heavy quark limit, in which the valence quark distributions 
become 5-functions centered at a; = 1/3, the moments {x'"-)^i-d approach 1/3". An 
extrapolation formula which explicitly satisfies both the heavy quark and chiral 
limits can be writteniEl 

{x")u-d - a« 1 + CLNAm^log „ ^ „ + bn ' (3) 

in which the coefficient of the leading nonanalytic (LNA) term, clna = ^(1 + 
3g^)/(47r/7r)^, is calculated from chiral perturbation theory.clia The mass ^ reflects 
the scale at which the Compton wavelength of the pion becomes comparable to the 
size of the hadron (without its pion cloud). Previous fitso to the latticejnoments of 
u~d suggest a value /i « 550 MeV, which is consistent with values foundj in analyses 
of nucleon static properties. The 6„ term in Eq. (^, 6„ = 1/3" — a„ (l — /i^CLNA), 
is included in order to provide a linear dependence on m^, and the mass scale A„ 
is set to be 5 GeV for all n. 

In Fig. 1 we show the best fit to the lattice dataBi for the n = 1 moment 
of u — d as a function of . Clearly, an extrapolation based on Eq. (0) provides a 
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much better fit to the lattice data and experiment than a hnear fit. Similar results 
are foundEJ for the n = 2 and n = 3 moments. Unfortunately, because all of the 
lattice data are in a region where the moments show little variation with m^, it is 
not possible to determine fi from the current data, and within the errors both the 
lattice data and the experimental values can be fitted with /z ranging from ~ 400 
to 700 MeV. Data at smaller quark masses are therefore crucial to constrain this 
parameter and guide an accurate extrapolation. 
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Fig. 1. The n = 1 moment of the u — d distribution, with a best fit to datatl'u using Eq. (|3|) with 
fj, = 550 MeV (sohd), and /i = 400 and 700 MeV (upper and lower dotted), compared with a 
Unear extrapolation (dashed). 



3. Quark Distributions from Lattice Moments 

Having established the appropriate way to extrapolate lattice data on PDF mo- 
ments, one can now ask much information can the existing data provide on the x 
dependence of the PDFs. The reconstruction of the complete x dependence in prin- 
ciple requires infinitely many moments, however, on the lattice, because of operator 
mixing for operators with spin > 5, all calculations have so far been restricted to 
n < 3. Nevertheless, as shown by Detmold et al.j3 considerable information on the 
shape of the valence Uy — dy distribution can already be inferred from just the lowest 
four moments. 

The most efficient way to reconstruct the PDF from a limited number of mo- 
ments is to determine the parameters of the PDF parameterization by directly fit- 
ting to the moments. For the isovector valence distribution accurate PDFs can be 
reconstructed from the lowest four moments, using the standard parameterization, 

xqy{x) EE x{q{x) - q{x)) = ax'^ {1 - x)\l + ey/x + jx) , (4) 

where /? can be related through Regge theory to the intercept of the p Regge 
trajectory, and S is given by perturbative QCD counting rules. 
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Note that some care must be taken when attempting to extract information on 
valence quark distributions from both the even and odd moments. While the even 
n moments of w — d correspond to the valence distribution, Uy — dy, the odd n 
moments correspond to the combination Uy — dy + 2u — 2d. The difference between 
these represents the famous violation of the Gottfried sum rule. Given sufficiently 
many moments oi u ~ d, one can reconstruct both the valence Uy — dy and u — d 
distributions by fitting the even and odd moments separately. However, at present 
there exists only a single data point for the even moments, n = 2 (the n = 
point corresponds to normalization), which makes it difficult to obtain accurate 
information. To minimize the error associated with reconstruction of the u — d 
distribution, we subtract the values of the phenomenological moments oi d — u 
from the calculated odd moments. The correction to the n = 1 moment is < 10%, 
while for n = 3 it is a fraction of a percent. 
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Fig. 2. Moments of the valence — distribution (scaled by 3" ) at the physical quark mass, 
extracted from fits to lattice data using a linear extrapolation (diamonds) and Eq. (0) (stars). 



The resulting fits to the moments of the valence Uy—dy distribution are displayed 
in Fig. 2 for both the linear and improved extrapolation, Eq. (||), with the shaded 
region around the latter corresponding to a 1 cr variation of the fit parameters 
from their optimal values. For clarity we plot 3" times the moments, so that the 
horizontal line at unity represents the heavy quark limit. 

The corresponding distributions x{uy — dy) are displayed in Fig. 3. Once again, 
the lightly shaded region represents a 1 a deviation from the central values, while 
the darker band illustrates the spread between global PDF fits. A comparison of 
the distribution reconstructed using the improved chiral extrapolation with the 
phenomenological distributions shows reasonably good agreement. On the other 
hand, the linear extrapolation gives a distribution (scaled by a factor 1/2 in the 
figure) which has a significantly higher peak, centered at x ~ 1/3, reminiscent of a 
heavy, constituent quark-like distribution. 
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Fig. 3. Physical valence x{u^ — d„) distribution, extracted using Eq. (solid), and a linear 
extrapolation, scaled by a factor 1/2 (dot— dashed) . 



4. Regge Behavior and the p Regge Trajectory 

According to Regge tlieory, tlie exponent f3 wltich governs the small-x beliavior of 
the distribution u„ — dy is related to the intercept (~ 1/2) of the isovector, C-odd, 
p Regge trajectory, and indeed the best fitEl to parameterizations of global data 
gives (3 ~ 0.48. The dependence of /? on the quark mass obtained from the fits to 
moments in Sec. 3 therefore allows one to predict the dependence of the Regge 
intercept. 

In addition to the intercept, one also needs to determine the slope of the tra- 
jectory as a function of niq. In the infinite mass limit, orbital excitations of mesons 
become energetically degenerate with the L = state. Within Regge theory, this 
is possible only if /? — s- oo as niq — > oo, which is consistent with the valence dis- 
tribution approaching a (5-function. One expects, therefore, that the slope should 
increase as rriq increases from its physical value. Using the lattice data for the p 
meson and the values of f3 generated from the fits in Fig. 2, one can then make 
predictions for the behavior of the masses of the orbital excitations as a function 
of rUq. 

In Fig. 4 we show the predicted p Regge trajectory at ttIt^ = 0.785 GeV, cor- 
responding to the strange quark mass, compared with the trajectory at the phys- 
ical light quark mass. A fit through the central values of P and the p massQ at 
ttItt = 0.785 GeV yields a slope which is larger than that of the trajectory at the 
physical quark mass, consistent with the expected trend towards the heavy quark 
limit. 

Although lattice data for the masses of orbital excitations are scarce, there have 
been some pioneering calculations of the 02 and pa meson masses by the UKQCD 
Collaboration,! indicated by the filled boxes in Fig. 4 (we use the fact that the a2 
trajectory lies on top of the p trajectory). Comparing with the predictions from 
the PDF analysis, the calculated p^ meson mass lies within the predicted band, 
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Fig. 4. Regge plot of the spin, J, versus t = (mass)-^ of mesons on the p trajectory, at the 
physical pion mass (dashed) and at m,r = 0.785 GeV (solid). The physical masses of the p, 02 and 
P3 mesons are indicated by stars, while the boxes represent lattice masses.Q The darker shaded 
region represents the statistical error in the extrapolation, while theJighter region indicates an 
estimate of the systematic error associated with the fitting procedure.^ 



albeit within large errors, while the 02 mass lies on the edge of the predicted range. 
Needless to say, further exploration of the masses of excited mesons within lattice 
QCD would be very helpful in testing these predictions. 



5. Helicity Distributions 

Lattice simulations of the spin-averaged quark distributions provide one of the 
benchmark calculations of hadron structure in lattice QCD. Having established 
confidence in the reliability of the lattice calculations through the correct treat- 
ment of the chiral behavior of the moments, we can now turn to other PDFs of 
the nucleon. Of particular importance are the helicity distributions, Ag(a;), which 
describe the distribution of the nucleon spin amongst its quark constituents. 

While a complete determination of the helicity distributions requires calculation 
of the singlet dis t ribution. E^Ag, for which there have only been exploratory lattice 
calculations,t2lllil a more basic challenge remains to understand the axial vector 
charge, 5^, given by the n = Oj^aoment of the isovector u—d distribution. The results 
of several lattice calculationsD'EHlil of gA are compiled in Fig:. 5. (Not included are 
results from recent simulations using domain wall fermionSjEj which appear to have 
strong finite volume dependence - see below.) When extrapolated linearly in vriq 
to the physical quark mass, the results are ~ 10-15% lower than the experimental 
value. Simulations with dynamical fermionsi (SESAM) yield results consistent with 



the earlier unquenched calculationajl 

On the other hand, the behavior of the moments (Aa;")„_ci is known in the 
chirall and heavy quark limits, and can be used to constrain the extrapolation 



Connecting Structure Functions on the Lattice with Phenomenology 7 



form, as in Eq. (||): 



ml 



(Ax")„_. = a„ ^1 + c,,,mi log j + 6„^^^ , (5) 

where clna = -(1 + 2.g^)/(47r/^)2, and 6„ = 5/3"+^ - a„(l - m^clna)- Using the 
same values of fi and A„ as in Fig. 1, the results of the chiral extrapolation using 
Eq. M are shown in Fig. 5. 



1.50 




Fig. 5. Lattice dat, 
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on the nucleon axial charge, gA, extrapolated linearly in (long 
: best fit (solid), and fits to extrema of error bars (short dashed). 



The result for gA is a downturn with decreasing m^, making the discrepancy 
with the empirical value larger. The reasons for this could be severalfold. Firstly, 
the chiral behavior in Eq. (J^ arises only from fluctuations N ttN N, whereas 
it is well known that the A plays an important role in gA through the fluctuation 
N ttA N. On the other hand, in the chiral limit the A has been shownB 
to give vanishing chiral contributions at order log , suggesting that higher 
order effects in the chiral expansion are likely to be important. The inclusion of 
these effects is currently under investigation. 

More importantly, perhaps, to fully incorporate the pion cloud in a lattice simu- 
lation a sufficiently large volume must be used. Indeed, there are some indicationsllHl 
of large finite volume effects for gA, which tend to increase gA in comparison with 
the results from smaller volumes. Clearly, in light of the results on the chiral ex- 
trapolation, it is imperative to perform simulations on larger lattices to understand 
the source of the discrepancy. 

6. Conclusion 

In this report we have highlighted the importance of model independent constraints 
from the chiral and heavy quark limits of QCD in the extrapolation of lattice data 
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on parton distribution moments. Inclusion of the nonanalytic structure associated 
with the infrared behavior of Goldstone boson loops leads to a resolution of a long- 
standing discrepancy between lattice data on low moments of the spin-averaged 
u — d distribution and experiment. 

The importance of ensuring the correct chiral behavior is further illustrated by 
comparing the x distributions obtained by extrapolating the lattice data using a 
linear and a chirally symmetric fit. While the latter gives an x distribution which is 
in quite good agreement with the phenomenological fits, the linearly extrapolated 
data give distributions with the wrong small- a; behavior, which translates into a 
much more pronounced peak at a; '--^ 1/3, reminiscent of a heavy, constituent quark- 
like distribution. Our analysis suggests an intriguing connection between the small- a; 
behavior of the valence distributions and the ruq dependence of meson masses on 
Regge trajectories, which should be tested more thoroughly in future simulations 
of the excited hadron spectrum. 

Finally, we have highlighted the need for further study of moments of the helicity 
distributions, and the axial vector charge of the nucleon in particular, which appears 
to be underestimated in lattice simulations. The inclusion of the pion cloud of the 
nucleon leads to a larger discrepancy at the physical quark mass, indicating that 
lattice artifacts such as finite volume effects may not yet be under control. Further 
data on larger lattices, and at smaller quark masses, will be necessary to resolve 
this issue. 
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